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Abstract 

It is shown that the series of renormalon-type graphs, which consist in 
the chain of insertions to one soft (hard) gluon(photon) hne is in fact ill defined. 
Each new type of insertions, which appears in the higher orders of perturbation 
theory, generates the correction to renormalon of the order of ~ 1. However, 
this series of the corrections to the asymptotics although have no small param- 
eter but hopefully is not the asymptotic one. The consideration based on the 
use of the renormalization group equation for effective charge is supported by 
the direct diagrammatic picture. 
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1. For many years it was generally accepted that the true asymptotics of the 
perturbation theory in theories with running coupling constant is determined by the 
renormalons [0, |^. Nowadays, the renewed interest is demonstrated in asymptotic 



estimates of perturbative series [3-11]. It results even in the recent attempts to 



use the renormalon for direct calculation of experimentally measurable quantities. 

However, the accurate determination of renormalon-type asymptotics appears to 
be not so simple problem. It was recognized by the speciahsts [|, ^ || that the overall 
normalization factor of the renormalon could not be found without taking into ac- 
count of all terms of the expansion of, say, the Gell-Mann-Low function. However, the 
usual proof of this fact do not refer on the direct counting of the Feynman graphs (see 
also the discussion after our equation ([T^)). Also recently Vainshtein and Zakharov 
||T0| have found the new source of uncertainty, which makes difficult the quantita- 
tive finding of the renormalon asymptotics (at least for the ultraviolet renormalon). 
The generally considered renormalon chain of graphs is formed by dressing of one 
gluon(photon) line by various insertions. The authors of ref. 0] have broken this 
tradition and considered the contribution to the ultraviolet renormalon of diagrams 
with two, three etc. dressed lines. Their analysis based on the operator product 
expansion, showed that contributions to the asymptotics from diagrams with any few 
dressed gluon lines are all of the same order of magnitude and also are parametrically 
larger than the usual renormalon with only one dressed line. 

In the present paper we would like to continue the analysis of the renormalon-type 
chain of Feynman diagrams. We consider the dressing of single gluon(photon) line as 
for the traditional renormalon, but try to estimate the role of the arbitrary high order 
insertions to this "dressed" gluon. It will be shown that each new type of insertions 
generates the correction to renormalon of the order of ~ 1. Thus in order to find 
the overall normalization of the asymptotics one has to calculate all the coefficients 
bo, bi, 62, ••• of the expansion of the renormalization group equation for, say, effective 
charge (see the eq. @ below). Our approach is equally valid for both infrared and 
ultraviolet renormalons. However, because there exists another source of problems 



for the ultraviolet renormalon [|I0|, we will concentrate our attention on the infrared 
one. 

It is a tradition now to consider the renormalon for QED. In this paper we 
will also discuss only the QED-type diagrams of the perturbation theory, with- 
out the self-interaction of gluons. However, while considering the renormalization 
group coefficients we would like to use the QCD values 60 = 4^ (y^c — fA'^/) , 

bi = (^j^y {y^c - f^/^c + Nf/Nc^ , ... . Thus by hand we avoid the problem 
of the gauge-independent definition of the renormalon chain in QCD. 

2. The contribution of the diagrams with exchange of one soft gluon(photon) to 
some "physical" quantity (say Re+e-^hadrons , or the correlator of two currents) has 
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the generic form 

R= f . (1) 

This is the first infrared renormalon, because of the contribution proportional to 
dk'^/Q^ vanishes due to the gauge invariance. The Feynman graphs corresponding 
to this value are shown in fig. 1. More precisely the graphs of fig. 1 correspond to 
(0) with the running coupling constant a{k) replaced by the fixed value ao = a{Q) ■ 
In we have written down the effective running coupling constant a{k) = a^ffik) 
instead of Oq. a^jfik) is trivially connected with the transverse part of the gluon 
propagator G^^{k). The function a{k) satisfies the renormalization group equation: 

^ = boa^ + bia^ + b2a^ + ... , (2) 
ax 

^ = ^"(^) ' 

where for practical applications we will use Nf = Nc = 3 , Anbo = 9 , IGn'^bi = 64 . 
Starting from 62 the coefficients of the expansion (H) depend on the renormalization 
scheme. It is to be noted here that we have fixed the renormalization scheme by 
considering the effective charge. Thus our coefficients 62,63,... are neither the free 
parameters, nor the known, say, for MS scheme, b2{MS), 63 (MS") . 

In order to find the leading contribution to the asymptotics, one may neglect 61, 62, 
etc. in (^. Now instead of 



9 roo roo p 



R = — / a(x)e dx = / 2dx 



-2x 



ao Jo Jo 1 — bottox 

N=0 N=0 \ ^ / 



N\ . (3) 



Here the first equality is also the explicit definition of the quantity, which we are 
going to consider in this paper. 

It is seen immediately from @ that if 60 > 0, our renormalon is ill defined due to 
the Landau pole. However, in this paper we will consider only the asymptotics of the 
perturbation theory and will not concern the issue of the nonperturbatve definition 
of the integral @. 

The way of calculation of the renormalon now seems straightforward. One should 
find, step by step, the function a{x) from @, substitute it into (H), and look for 
the asymptotics. Nevertheless, before passing to this formal way, let us illustrate 
the role of complicated contributions to the renormalon by the explicit estimate of 
corresponding Feynman graphs. 
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The fig. 2 shows the chain of diagrams corresponding to renormahzation of the soft 
gluon fine in the leading approximation (|^). As we have said above, we show only the 
QED-type diagrams consisting of fermionic bubbles without gluon self-interaction. 
In the A^-th order of perturbation theory each of the bubbles from fig. 2 generates 
the factor 

6o«oln(^) (4) 

in the integrand of Now the difference between theories (QCD,QED) is hidden 

in the factor bo, accompanying the single bubble. 

Now let us replace two of the simple bubbles by the more complicated two loop 
diagram, as it is shown in fig. 3. Again we have shown only the QED diagram. 
Moreover we have shown explicitly only one of the two diagrams (compare with 
fig. 1) of the second order in oq. As it is indicated on the figure, the - two loop 
bubble generates the factor 

M^ln(^) (5) 

in the integrand, which has one power of large logarithm less (or one cto more) than 
the leading order contribution (^). On the other hand a large combinatorical factor 
A^ — 1 appears due to a number of permutations of the second order bubble among 
the simple bubbles. As a result one has 




(6) 



Thus we can see that taking into account one second order insertion into the soft 
gluon line leads to the correction of the order of one to the trivial asymptotics (|^). 
Summation over the number of the second order insertions shown in fig. 3 leads to a 
simple exponentiation of this correction 



In last few years it became very popular [Tl[| to consider the renormalon in 
the limit Nf — > oo. From this point of view the correction (|^) is nothing more than 
the ~ 1/Nf correction, because of 26i/6q ~ 1/Nf ^ 1 . However for practically 

L28 906 678 
81 ' 625 ' 529 ' 



interesting A'c = 3 and Nf = 3, 4, 5 one has, respectively, 26i/6g - ^'^^ 



Consider now the more complicated diagram of fig. 4 with dressing of the internal 
gluon line of the second order bubble. At this point it is natural to write down 
explicitly the last integration over internal momentum of the two loop diagram 

'Q'\ ,2 f^'dq^ 
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Now the dressing of the gluon hne evidently leads to 



hial 



6o«oln 



g2 



n + 



6o«oln[^ 



(9) 



Thus up to the overall factor the contribution of diagram of fig. 4 coincides with 
that of fig. 3. Summation over the number of simple bubbles inserted into the large 
bubble - n naturally leads to In(A^). Also taking into account a number of large 
bubbles of fig. 4 allows to exponentiate the correction 



N 



A^lexp 



'boao 



N 



(10) 



This is the generally recognized expression for the infrared renormalon. Our argu- 
ments up to this stage simply repeat the line of reasoning of the paper though 
may be in more details. 

It is clear, that the argument of the exponent in (^) was found with the ~ 1/ In(A^) 
accuracy and therefore the nontrivial overall factor as well as the function of A^, weaker 
than N'^, may appear in ([l^), as we consider in detail in the following section. 

Now let us consider the three loop correction to the renormalon, shown in fig. 5. 
Like it was done before, we show only one example of the third order diagram. All 
these contributions generate the factor in the integrand of (^ 



Thus here we have two extra ao, which at first glance could not be compensated by 
one combinatorical factor A^ and hence the diagram of fig. 5 seems to generate only 
the ~ 1/A^ correction to renormalon. In particular such conclusion was drown by 
Zakharov 0. However, let us see, what happens if one dresses the internal gluon 
lines of the three loop diagram (fig. 6). Now summation over the number of trivial 
corrections inserted into a large diagram ni,n2 gives: 



6200 In 



62^0 In 




n (N-rii 



n2 



X (aoA^) ~ &2«o In 



(12) 



Here the factor [rii + 722 + 1)^^ in the sum appears after integration over the internal 
momentum of the large bubble (see (^), while the factor (A^ — ni — n2 — 2) accounts 
for the number of permutations of our large bubble with simple small bubbles on the 
main gluon line. So we see that after dressing of all gluon lines the three loop (~ 62) 
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diagram generates the correction to renormalon of the order of ~ 1. One can easily 
show that four loop (~ 63), five loop (~ 64) etc. diagrams generate the corrections 
of the same order of magnitude. Previously the analogous proof of the importance of 
the high loop corrections to the renormalon chain was done by Mueller [0, but this 
result was not published. 

3. In this section we would like to turn back to the direct analysis of the renor- 
malon (|1]J^) without referring to the Feynman diagrams. It is easy to integrate for- 
mally the renormalization group equation 



1 



1 



a \ 



\ 7- In — - C2(a - ao) - C3- 

a ao bo KaQ/ 



a 



an 



where, Cs = - bj/bl , C3 = 63/60 - 26261/60 + ^i/^l 
introduce the new set of variables 



t 

boa 



2x 



a 



6o«o 
2 

26i 
' IF ' 



- ... = 60a; , (13) 
Also it is convenient to 

(14) 



/3i = /3 

Now the renormalon (^) takes the form 
R - ' 



— e at , 
a 



(15) 



a 
a 



1 — at — /3ia\n{a/a) — [32{a — a)a — I3z{ci^ — a^) 



a^ja — 



Although we have shown, that the corrections to renormalon generated by the 
high order contributions to the renormalization group equation 62^'^, 630^... (0) are 
not small, the correction induced by the second term bia^ still plays an outstanding 
role due to the additional enhancement by In(A^) (see (p!oD). Therefore one should be 
more careful, while treating these enhanced contributions. Moreover, the iterations 
of the second order diagram of fig. 3, e.g. as it is shown in fig. 7, may also have some 
additional enhancement. Thus at the first stage let us omit ^2,^3, Pa ... in ( p!5[ ) and 
consider the truncated effective charge: 



a 
a 



1 



1 — at 



pa m(a/aj 

This is the transcendental equation for the function a = a{a), which may be solved 
iteratively. For estimate of the A^-th order of perturbation theory we will often use 
the formula for A^-th term of the expansion of the integral in powers of a 



'dt 



1-at 



Pa 



1 — at 



In 



1 - at 



kl 



o 



InN 



(17) 



N 
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Here both m and k are supposed to be large m,k In(iV) . Everywhere in this 
section we suppose that not only the A^, but also the In(iV) is a large parameter and 
neglect all the corrections of the order of ~ l/ln(A^). In order to derive the formula 
(p!7|) one has to use the asymptotics of gamma-function as well as the trivial identity 



Mp)T 



hm (^^ 



p 



Now it is easy to make the first iteration in (|rBl|TH) 



{Ri} 



N 



1 — at + /5aln(l — at) 



1 - at 



E 



Pa 



In 



1 — at 1 — at 



(19) 



1 



p 



'6o«o 



N r 



N 



n/3 



plnN 



N\ 1 + 



InN 



Here for brevity we "forget" to indicate, that only the A^-th term of the expansion of 
both integrals over a is considered. The asymptotics (|19D should be compared with 
(10). One can see that the consistent treatment of the diagrams of fig. 4, which in 
fact we have done in (|l^), results in the nontrivial small factor (InA^)"'^ as compared 
to the naive result (p^O]). 

Consider now the effect of iteration of two loop correction (fig. 7), which in terms 
of a means 



^ Oi / second iteration 

To this end one has simply to replace in (|T9| 



1 - at + Pa\n{l - at + ^aln(l - at)) 



(20) 



Pa 



1 — at 

Pa 



In- 



1 



at 



-Pa\n{l - at + /51n(l - at)) 



at 



L 



In 1 - 



Pa 



1 — at 



at 



(21) 



where we introduce L = In (jr^)- As we have seen from ( |T^JT^ ), from the point 

of view of perturbation theory ~ (InA^)^^, while L ^ InA^. As a result the 
expression in curly brackets in ( pT] ) is of the form 



L 



In 1 



Pa 



at 



-L 



1-0 



'ln(lnAr))" 
InA^ 



and may be written as the exponent of a small quantity. Thus instead of 
gets 





p 


\i l] 


f 
L 

= E 

n=0 




1 - at 




1-at 


1 — at 



p+n Z_i 



T{p/L) 



(22) 



one 



(23) 
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In terms of Feynman graphs the variables p and n are respectively the number of 
large second order bubbles and number of internal second order bubbles (see fig. 7). 
Also all wavy lines in fig. 7 are supposed to be renormalized by arbitrary number of 
simple one loop insertions (like in fig. 2). 

Substitution of (|2^) into (|19|) gives 



q=On=0 y- 



/31n 



q 3^-1 



^ \lnN J 



(24) 



Here we have denoted p + n (p3D by q. Due to [pinNY/q\ the series in q ( p4D has a 
narrow peak. Therefore, up to ~ 1/ In iV corrections 



N 



where x = n/ In(A^) 
In(lnA^) > 1 





'N' 


^ f 


)■ 




N\ / 




J. 


Jo 



^l3-x-l^-x\n{lnN) 



dx 



T{(3-x) 



(25) 



). The explicit integration here may be performed only if the 



{R2} 



'6otto 



N 



N 



N 



l3\n{\nN) 



/3 



mi 1 + 



In(lniV) 



(26) 



Thus again we can see (compare (|26| ) and ([l9|)) that taking into account the diagrams 
of fig. 7, or in other words, making the second iteration (PD| ) in the transcendental 
equation ([T6|), leads to parametrically large renormalization of the renormalon. 

In the analogous way one may perform the third iteration of (0). The result 
reads 



'&o«o 



N r 



TV 



N 
J 



m / dxi 



dxo 



(xi -xa)^-"^-^ x^2'~"'~^ 
T{(3-xi) r(a;i-a;2) 



-X2 In(lnAf) 



(27) 

Here the integration may be performed analytically only if the ln(ln(lnA^)) ^ 1 
resulting in the same expression as (^), but with In(lnA^) replaced by ln(ln(lnA^)) 
both in the main formula and in the estimate of the error. 



For the case of arbitrary k + 1 iterations in (^) the generalization of 
(p7|) leads to 



and 



{Rk+l} 



N 



N 



iXi — X2Y 



N 
J 

[X2 - X3) 



n/3 



A^! / dxi / dx2 ■■ 



Xk-l 



dxu 



(3 -xi 



T{l + (3-xi) 



(28) 



[Xk-l - Xk 



\Xk-2—Xk-\ 



'^fc — l Xk 
Xi, 



-Xk In(lnAf) 



r(i + xi - X2) r(i + X2- xs) "' r(i + Xk-i - x,] 
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Again one can calculate this integral analytically if the A;-th logarithm is large 
ln(ln( ... In A^) ... ) S> 1, though the expression (28) itself has the accuracy ~ 1/ In A^. 



Looking at the formula ([28|), one may even doubt, whether it has any finite limit 
for large k. In order to prove this finiteness consider in more details the last integral 
in (|28|) 



" rfXfc ^^;^ g-x.ln{ln7V)^ (29) 

Vil + Xk-i-Xk) Xk 



Xfe_2-a;'fc_i-l / -I , 
Xk-1 



1+ CmnlXk-2Xk-l{^'^Xk-iy] . 



m+n>0 

Substitution of this result into (|28D naturally leads to the following estimate 

const^^ ^ 



{i?fc+i}^ = {i?fcV(l + 0(— 1 . (30) 



By the way we have shown here that the integral in ( |28D at large k forgets effectively 
about the value of In(A^) . 

Thus we arrive at the surprising result. If one performs only one (a few) iter- 
ations while solving the transcendental equation for the effective charge ([16|), the 
asymptotics of the perturbation theory differs drastically (|19,25-27) from the gener- 



ally recognized renormalon (pDI) . Nevertheless, after taking into account the infinite 
number of iterations (^) the renormalon ([T0|) is restored up to some overall constant. 

4. In the present paper we have no plan to perform the complete analysis of the 
contribution to renormalon from the high order terms of the renormalization group 
equation @) h2a'^,h^ct' ... (in particular because the coefficients 62,^3 themselves are 
not known). It will be enough for us to show that all the corrections to renormalon 
due to these high order terms are equally important. To this end we will consider the 
contribution to renormalon from the multi-loop corrections in the equation (^ only 
in the leading nontrivial approximation. 



As we have seen (p!^-[T5|) in the formal solution of the renormalization group equa- 
tion the coefficients 6„ have been converted into = 2"6„/6o^^ + .... In the linear 
in ^2,^3, ■■■■ approximation the effective charge ([T5|) takes the form 



a «! 
a a 



1 + ^2^1 + (33^1+ ...) , (31) 



where ai is the solution of the equation (|T6D. While going from (|15|) to (^T]) we have 
ignored the /32a^ in comparison with /?2«^, (^^a^ in comparison with Psa^ and so on. 
With the experience of the previous section one can easily show, that the omitted 
contributions are of the order of ~ 1/A^. 
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Now it is easy to repeat all the logics, which led us to the equation 



{Rk+i} 



N 



N 



N' 
J 



/3 



/3 



A^! / dxi (32) 







j^l + _ + - x^f + ... ^ £ dx, ... , 

where the part of the formula after integration over dx2 simply repeats the corre- 
sponding part of the equation (^) . Keeping in mind the finiteness of both (pSD and 
(^) (see ( p9D and (|30|)) one can see that all /92, /?3, ... make contributions of the order 
of ~ 1 to the renormalon. 

Nevertheless the formula (^2]) allows one to draw not only the pessimistic con- 
clusions. All the high order (~ terms in ( P^ ) have appeared in the combination 
I3n/n\ . Thus one may hope, that at least the series of the corrections to the renor- 
malon is not the asymptotic one. 

5. The formulas like (|28|) still are too complicated for practical use. In this section 
we would like to develop another method for calculation of the renormalon-type 
asymptotics, which at least enables to perform the simple numerical computations. 

Let us again introduce the new variables instead of (|H|) 

Here in the last equality we use, that in all corrections to the trivial renormalon chain 
(|) the variable t may be replaced by the saddle-point value t = up to corrections 
of the order of ~ [lnA^]^/A^. Now the truncated formula for the effective charge (16) 
takes the form 

/ = ^—r^ , (34) 

where 7 ~ 1/A^ . The A^-th term of the expansion of / in the series over s just gives 
the asymptotics 

f = l + ^C„(7)s", (35) 

1 

{R}n = mc^ip/N) . 

It is easy to show that f{s) satisfies the differential equation 

r! df /(/ — I) , 

J = -r = + ^^JJ ' (36) 

as s 
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which leads to the recursion relation for the coefficients of the expansion 

^ n=l ^ n=l 



(37) 



with the initial condition Ci = 1. For analytical investigation this formula is even 
more complicated, than the original equations ([T5|), (p!6|). On the other hand, it is 
easy to solve (|37|) numerically. 



We have performed the numerical simulations with 60 = 9, 61 = 64 and N up to 
16000. The result reads 



i? = 0.275](^^ao 



N 



(3^ 



Let us recall that this result is obtained for = /^s = .. = (see (p!^),([T5|)) or, in 
other words, we were interested in the solution of the equation (compare with (0)) 



da 
dx 



hi 
bo 



(39) 



a 



If one expands the r.h.s. of this equation in series in a, the terms proportional to 
b1,bl, ... will mimic effectively the ~ 62, ^3, ••• terms of the full equation (0). On the 
other hand, one may consider the pure truncated renormalization group equation: 



da 
dx 



In this case the formulas 

T 

1 — s — 7sln 



/ 



boa'^ + bia^ . 
|36|) are modified 



(40) 



(41) 



More precisely this formula for the effective charge corresponds to the approximate 
solution of (1^), but the omitted terms lead only to the ~ 1/A^ corrections to renor- 
malon. Again it is useful to expand f{s) in the series (0). The corresponding 
recursion relation takes the more complicated form 



C, 



k+l 



1 " 

T / , CnCk+l- 



(42) 



k 



Ck + 2 ^ CnCk-n + ^ CnCmCk-n-m 

~n—m—l 
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where again Ci = 1 and all C-s with nonpositive numbers are supposed to be zero. 
After numerical simulation of about 1000 terms of the series (42) we have found 



i? = 0.13^f^«oj iV^iV! . (43) 



All the difference between the two definitions of the effective charge (^9|) and ( ^0|) 



lies in the different choice of higher order corrections to the renormalization group 
equation 62,^3 ••• • Thus the two asymptotics (|38|) , (|43|) show exphcitly that these 
corrections make the contribution of the order of ~ 1 to the overall normalization of 
renormalon, in complete agreement with the conclusion of the preceding section. 
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Figure captions 



Fig 1. The renormalon-type graphs with exchange of one soft gluon. The 
internal gluon line will be dressed in the following figures. 

Fig 2. The example of the simplest chain of diagrams corresponding 
to the renormalization of the soft gluon line. Each bubble generates the factor 

Mo In (gv^')- 

Fig 3. The example of the chain of diagrams with one two - loop insertion in 
soft gluon line. This second order bubble generates the factor bia^hi (^Q^/k'^^. 

Fig 4. The dressing of the internal gluon line of the second order bubble, 
shown in fig. 5, by n simple bubbles. 

Fig 5. Three - loop insertion to soft gluon line. This diagram generate the 
factor &2Q;oln (^Q^/k'^^ and makes only 1/N correction to the renormalon. 

Fig 6. The same as in fig. 5, but with dressing of two internal gluon lines by 
the simple chains of bubbles. The summation over rii and n2 allows to compensate 
the extra as of the contribution of fig. 5, thus leading to the correction of the order 
of ~ 1 to the renormalon. 

Fig 7. An example of the graph which corresponds to the second iteration 
of two - loop correction (see fig. 3). Here p is the number of (large) dressed bubbles 
which lie along the soft gluon line and n is the total number of internal (small) dressed 
bubbles. 
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Fig.l 



bo aQ ln(Q /k ) 




Fig.2 

■ ln(Q7k ) 




Fig.3 



Fig.6 




Fig.7 



